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We study the finite size efi^ect of the Kondo screening cloud in a double-quantum-dot setup via 
a large-N slave-boson mean-field theory. In this setup, one of the dots is embedded in a close 
, metallic ring with a finite size L and the other dot is side-coupled to the embedded dot via an anti- 

■ ferromagnetic spin-spin exchange coupling with the strength K. The antiferromagnetic coupling 

' favors the local spin-singlet and suppresses the Kondo screening. The effective Kondo temperature 

Tfe (proptotional to the inverse of the Kondo screening cloud size) shows the Kosterlitz-Thouless 
(KT) scaling at finite sizes, indicating the quantum transition of the KT type between the Kondo 
screened phase for K < Kc and the local spin-singlet phase for K > Kc in the thermodynamic limit 
with Kc being the critical value. The mean-field phase diagram as a function of l/L and K shows 
a crossover between Kondo and local spin-singlet ground states for K < Kc {L — 4n, 4n -I- 1, 4n -I- 3) 
and for K > Kc {L = An + 2). To look into the crossover region more closely, the local density of 
states on the quantum dot and the persistent current at finite sizes with different values of K are 
also calculated. 
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I. Introduction 



^3 ; Kondo effect [H, the screening of magnetic impurities by the conduction electron reservoir, has been intensively 
y ■ studied over the last half-century. Recently, there has been revival interest in Kondo effect in semiconductor quantum 
^ ' dot devices due to the progress in fabracating nanostructures 3j. When there are odd number of electrons on the 
^ ' dot, Kondo effect overcomes the Coulomb Blockade, leading to a narrow pronounced resonance peak in local impurity 
Q density of states and enhanced conductance through the quantum dot[2|0]- The Kondo effect is characterized by the 
O single energy scale , the Kondo temperature Tk ~ De~^^"'^, below which the Kondo screening develops. Here, D, 
J in Tk are the conduction electron bandwidth and the dimensionless Kondo coupling, respectively. In a finite-sized 
psj . mesoscopic device, Tk sets the length scale ~ Hvp/Tk associated with the size of the Kondo screening cloud where 
^ ' the cloud of electrons with the size of order of surrounding the magnetic impurity form spin-singlet state with 
On : it0]. Here, vp is the Fermi velocity. For typical values of Tk, is around 0.1 l/im, which is comparable to the 
■ typical size of the quantum dot device, leading to the finite size effect of the Kondo screening cloud. This effect 
. has been investigated in a quantum dot embedded in a mesoscopic ring threaded by a magnetic field [4i][5j^]. The 
experimentally measurable persistent current induced by the magnetic flux has been shown to be sensitive to the ratio 
of the size of the ring L and C° HiS- As L < the Kondo screening cloud does not develop completely, giving 
' rise to the suppression of Kondo effect and a reduction of the persistent current even if the temperature is lower than 
Tk; while as for L ^ the Kondo screening cloud is formed, leading to enhanced persistent current. By measuring 
^ the transport properties (such as: persistent current) as a function of the system size, we gain insights on how the 
Kondo screening cloud is formed as the system size is increased. This idea offers another route to realize Wilson's 
' Numerical Renormalization Group (NRG)[7| idea on the Kondo problem. 
^ . Very recently, study of the Kondo effect has been extended to the coupled double-quantum-dot setups via antifer- 
■ ■ ' romagnetic RKKY interactions [8] [9] jl^l [HI- The Kondo effect in such double-dot systems competes with the RKKY 
interactions, giving rise to quantum phase transition in the context of the well-known two-impurity Kondo problem 
between the Kondo and local spin-singlet ground states [l^. Close to the quantum critical point physical observables 
at finite temperatures exhibit non-Fermi liquid behaviors. The crossover behaviors between these two phases can be 
accessed by changing temperatures. However, at zero temperature but at finite-sizes, the size of the Kondo screening 
cloud provides us with an alternative route to the crossover behaviors in the two-impurity Kondo problem[l^. 

In this paper, we study the Kondo effect in a side-coupled double quantum dot system embedded in a finite-sized 
mesoscopic ring. A similar side-coupled double-quantum-dot system has been studied where one of the dot is coupled 
to Fermi-liquid leads of conduction electrons with continuous spectrum[l3]. Unlike the two- impurity Kondo system, 
in the side-coupled double-dot system the Kondo phase is fragile and unstable towards the local spin-singlet state for 
any infinitesmall antiferromagnetic spin-exchange coupling, and the transition is of the Kosterliz-Thouless type. Note 
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that a related setup consisting of double quantum dots in parallel connected to a mesoscopic ring shows a quantum 
phase transition of the pseudogap Anderson model when the magnetic flux is tuned |15l|. Here, we are interested in not 
only the nature of the Kondo-to-spin-singlet quantum phase transition in our proposed setup but also the crossover 
phenomena via a systematic study of the system at finite sizes. The inverse of the system size 1/L plays a similar 
role as temperature T . The finite temperature crossover behaviors between the above two quantum ground states can 
therefore be accessed effectively via changing the size of the ring. 

We focus on the following three quantities to investigate this issue: the effective Kondo temperature Tk (inversely 
proportional to the size of the Kondo screening cloud in the presence of the antiferromagnetic spin exchange 
couphng), the local density of states (DOS) pqd{i-o) on the quantum dot, and the persistent current (PC) /. The plan 
of the paper is as follows. In Section II., we introduce the model and present a large-N mean-field treatment of the 
model. In Section III., we present our results on the effective Kondo temperature, the local DOS on the dot and the 
persistent current. We also give detail explanations of our results. The conclusions are given in Section IV.. 




FIG. 1: Sketch of the model. tL/ta is the hopping coefEcient between dot 1 and its left/right conducting island on the ring, 
is the magnetic flux (in unit of 27r/$o with <E>o = hc/e) through the ring, K is antiferromagnetic coupling between dot 1 and 
dot 2. 



II. Model 

Our model describes a double-quantum-dot system embedded in electronic reservoir in a finite-sized ring (see Fig.). 
In this setup, only one of the dot (dot 1) is coupled to the ring via hoping; while the other dot (dot 2) is coupled 
only to dot 1 via antiferromagnetic spin-exchange interaction with a coupling strength K and is decoupled from the 
reservoir of the ring. The ring serves as electronic reservoir consists of L — 1 conducting islands. We consider here 
the half-filled case {Ne = L + \ where N,. is the total number of electrons on the ring and on the two dots) . The key 
point in this setup with a closed geometry is that the Kondo screening cloud is trapped in the ring and can not escape 
into the external leadsf^. Nevertheless, one can measure the transmission probability through the quantum dot 1 
embedded in the ring by measuring the persistent current [lit (^^^ details in Sec. III.). Here we assume a direct 
antiferromagnetic spin exchange coupling between two quantum dots which competes with the Kondo effect. Note 
that though it has been known experimentally that the antiferromagnetic spin-spin coupling (RKKY interaction) 
between a double-quantum-dot system can be induced naturally by a conduction electron island in the middle of the 
two dotsQ, the direct antiferromagnetic spin exchange coupling may in principle be generated via the second-order 
hoping process directly between the two dots if the two quantum dots are close enough to each other. Since in both 
cases we expect suppression of the Kondo effect due to local spin-spin exchange coupling, to simplify our calculations 
we assume a direct antiferromagnetic spin-exchange coupling between the two dots. 

The Hamiltonian of our model is given by: 



H ^Ha + KSi ■ S2 



(1) 
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where Ha represents a quantum dot embedded in a ring by the Anderson impurity model. 

L-l 

Ha ^ -tY. Et^^^+i + ^-^J + E + E (2) 

j = l <T i=l,2 (T,i=l,2 

(7 

where Cj, c?i, ^2 represent the annihilation operators of electron on site i of the ring, dot 1 and dot 2, respectively, 
the phase factor is defined by = 27r$/$o with $ being the magnetic flux going through the ring and $o = hcje 
being the flux quantum, t is the electron hoping within the tight-binding ring, and i^/ij are the hopings between 
dot 1 and the two neighboring ring electrons (ci_i/ci). Here, 5*1 and 5*2 represent spin operators of dot 1 and dot 
2, respectively with Si = X^cro-' d\^(Ja(T'di(,' ■ L is the total number of electrons in the system excluding dot 2. The 
well-known Kondo limit is reached for U ^ t. Here, we consider a simple limit in the Kondo regime, U oo where 
the slave-boson mean- field (SBMF) approached is applicable to simplify the quartic U term. 

Further progress can be made by decoupling the quartic spin-spin interaction 5*1 ■ S2 into qudratic one via the 
Hubbard-Stratonovish transformation in the framework of the large-N SU{N) mean-field theory where the symmetry 
of the Hamiltonian is generalized from the SU{2) with spin degeneracy being two to SU{N) with spin degeneracy 
being — > 00 1]. In the SU(N) generalization of the slave-boson representation, we have dia = fiab\, ct is the flavor 
of the spin: a = 1, 2, ...N, i = 1, 2 is the indix for the two quantum dots. The local constraints to enforce the single 
occupancy on dot 1 and 2 are given by 

Y.ftf^o.+blb,^^ (3) 

a 

The mean-field Hamiltonian for Ha is therefore given by, 

L-l N N N 

NHaMF = -tEE[4^J + l"+^-^-]+E^~rfA^"/la + E^^^/2a/2a (4) 
j—1 a— 1 Q — 1 a— 1 

N 

- Y.^tLfLcio.^tBeM^(^yl^iJio. + H.c] + \{hl - 1) + \{F - 1) 

a=l 

where ^L/i? = boti^/R with 69 being the expectation value of the bi boson on dot 1, 69 =< 61 >, b ~< 62 >, A 
(A) is the Lagrange multiplier which enforces the local contraint on dot 1 (2), id = — A, = — A. The quartic 



antiferromagnetic spin-spin interaction can be decoupled via the mean-field variable x 12|: 



X = — < d[ d2,a > (5) 



. Therefore, 
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if^i • ^2 = - E f /^/2" - H-c. + ^ (6) 

a = l 



The full mean-field Hamiltonian is given by 



N 2 

Hmf = Hamf - E ^/^/2« - H.c. + (7) 



a=l 

The mean- field energy Emf can be obtained by diagonalizing the above Hamiltonian and can be expressed as 



Emf = E ^™--(^' ^0) + ^ + Kbl - 1) + W - 1), (8) 
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where e^.o- are the eigenvalues of the Hamiltonian matrix in Hmf, the summation over m includes all occupied levels 
of Hmf- The values of the mean-field variables A, bo and x are determined by minimizing Emf with respect to bo 
and X- 

QEmf _ OEmf _ QEmf _ „ 

subject to the constraint ^^qI'^ — ^^gj^'' = 0. The ground state energy Egg corresponds to the global minimum of 
Emf{^, a, bo, b, x) which satisfies the mean-field equations. Note that the advantages of taking the large- iV mean-field 
approach are: (i) in the large- iV limit the solutions from the mean-field equations are exact though the physical system 
corresponds to A'' = 2, and (ii) at finite N, a systematic correction to the mean-field results is possible though it 
is beyond the scope of this paper. The mean-field phase diagram can be mapped out via the solutions of the above 
mean- field equations. 
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FIG. 2: Schematic mean-field phase diagram of the model for (a) L = 4n, An + l^An + 3 (case (I)) and (b) L = 4n + 2 (case 
(II)). The red line indicates the crossover with K being fixed and the green line indicates the crossover with L being fixed. 



III. Results 

Three physical observables are identified to investigate the Kondo screening effect of our setup at finite sizes: (i). 
Effective Kondo temperature {T^) , (ii). Local density of states on dot 1 (LDOS) (pQo{i^)), (iii). Persistent current(PC). 
Details are shown below. 

Before we present our new results, it is useful to summarize the behavior of the model at K = (without antifer- 
romagnetic spin-spin coupling), which has been intensively studied [3| 00]. The Kondo resonance of a quantum dot 
embedded in a mesoscopic ring strongly depends on the finite size L (mod 4) , and the magnetic flux $ threading the 
ring. In particular, both the Kondo temperature Tk and the persistent current / follow universal scaling functions 
of £,1/ L at a fixed magnetic flux. In the Kondo regime, all physical observables, such as: Kondo temperature T^, 
persistent current / are enhanced as the size L increases, but with different crossover bahaviors in all four cases of 
L — An, L — An + 1, L = An + 2, and L = An + 3. The magnetic flux dependence of persistent current I exhibits 
a symmetry between size L and L + 2: Il{4>) — lL+2{'t' + tt), indicating that adding the magnetic flux of $0/2 (or 
adding a phase tt) is equivalent to switch the behavior of the PC from a system with size L to L + 2. 

With the previous results in mind, we may discuss the general properties for K > 0. Due to the antiferromagnetic 
spin-spin couping, we expect in this case the competition between the Kondo and the local spin-singlet ground states, 
leading to the quantum phase transition. In fact, quantum phase transitions in double-quantum-dot systems with 
antiferromagnetic RKKY interactions have been intensively studied in recent years in the framework of two- 



impurity Kondo problem 13[ where the quantum dots couple to the conduction electron Fermi sea with continuous 
spectrum. Two types of quantum phase transitions have been identified in these systems: the phase transition with 
a quantum critical point (QCP) and the one of the Kosterlitz-Thouless (KT) type. The characteristic behavior near 
QCP is the observables power-law dependence on the coupling strength relative to the critical point; while for the 
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KT transition the crossover energy scale exponentially depend on the distance to the critical point. The former 
(QCP) type of the quantum phase transition is realized in the double-dot systems where each of the dot couples to an 
independent conduction electron reservoirQQ, and the critical point separating the Kondo from the local spin-singlet 
phase is the well-known two-impurity Kondo fixed point [l^. The latter (KT) type exsits in a side-coupled double-dot 
system where only one of the dots coupled to the electron reservoir [l^. The Kondo resonance becomes more fragile 
in the side-coupled system so that an infinitesmall antiferromagnetic coupling is sufficient to suppress the Kondo 
effect and leads to the spin-singlet ground state. We expect the similar Kosterlitz-Thouless transition to occur in our 
side-coupled double-dot system embedded in a ring. However, the mesoscopic ring in our setup consists of a finite 
number of tight-binding electrons (instead of a Fermi sea with continuous spectrum), the details of the transition 
might be different from those in Ref.[3| (see below). 



A. Mean-field phase diagram 

After solving the mean-field equations, we summarize our main results in the schematic mean-field phase diagram 
as shown in Fig. [51 We find indeed the transition between the Kondo and spin-singlet phases is of the Kosterlitz- 
Thouless type (see below). However, we find that the critical point separatiiig the two phases is not at zero as 
shown in the similar side-coupled double-dot system studied previously in Ref. [l4j but at a finite value: Kc > 0. 

There are three regions in the phase diagram, corresponding to different mean-field solutions: 

1. For small K we find fog 7^ 0, A 7^ 0, x = 0. In the thermodynamic limit this is the Kondo phase studied in 
Ref.iii. 

2. For large K, we find fog = 0, A 7^ 0. The ground state for L — > cx) is the local spin-singlet phase where 
antiferromagnetic spin-spin coupling completely suppresses the Kondo effect. 

3. In the intermediate values of K, we find 60 = 0, A 7^ 0, x 7^ 0. This corresponds to the crossover region between 
Kondo and spin-singlet phases indicated in the shaded region in Figl2l This region is defined cither by Kd < K < Kc2 
for a fixed size L where Kd and Kc2 are the boundaries between the crossover region and the two stable phases or 
by Lci < L < Lc2 for a fixed K where Ld and Lc2 are defined in a similar way as Kd and Kc2- 

Note that in all the above three cases, we find & = 0, A = e^. We then systematically investigate the crossover 
behaviors along the following two different crossover paths, depending on the finite size L (mod 4): 

Case (I) (Figl2] (a)) holds for L = An, 4n -I- 1 , 4n -I- 3 where we find the crossover occurs for K < Kc- 

Case (II) (Figl2] (b)) occurs for L = 4n + 2 where we find the crossover ranges exist mainly for K > Kc- 

In both cases, we find the crossover energy scale T* follows the behavior of the typical Kosterlitz-Thouless transition: 

T* - cfk exp [~7:fk/{K - K^)] (10) 

where = c'T^, and both c and c' are non- universal constants. 

We investigate the Kondo effect in our setup at finite sizes by either changing the size L at a fixed K (or following 
the red (vertical) line in Figl2|) or changing K at a fixed size L (or following the green (horizontal) line in Fig. [2]). 
The finite size scaling of Tk indicates that in the thermodynamic limit L 00, Kd and Kc2 converge to a single 
critical value K^. Kd = Kc2 Kc (see below). 



B. The effective Kondo temperature 

In the single Kondo dot system embedded in a 2D electron gas (2DEG), it is well-known that the physical properties 
follow universal functions of {T/T^) in Kondo regime where is the Kondo temperature of the single dot in the 
thermodynamic limit In a single quantum dot embedded in a mesoscopic ring with a finite size, the effective Kondo 
temperature can be defined as[l| 



(11) 
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, where is the energy of the dot 1, Egs and E^^ correspond to the energy of the full system and of the tight-binding 
ring with the open boundary condition, respectively. The effective Kondo temperature Tk is the energy gain due to 
the coupling between quantum dot and the ring, which corresponds to the energy associated with the Kondo coupling. 
In the thermodynamic limit {L — > oo), approaches to T°. For the relevant parameters in the K = Anderson 
model of the single dot-ring setup: t = 1, = = 0.4t, = -0.8t, we find « 0.0189t, ~ 106. Note that for 
simplicity, we set t = 1 throughout the paper as our unit. 

In the presence of the antiferromagnetic spin-spin coupling the effective Kondo temperature Tk at a finite size in 
the context of our large-N slave-boson mean-field approach can be generalized from Eq. [TTlto: 



Tk — Cd - {E, 



e: 



gsl 



X2 

K 



(12) 



Note that the effective Tk defined in EqUT] and Eq[l2] are slightly different from that defined in Ref . [3] where Tk is 
measured relative to the hightest occupied level defined as ef, corresponding to the Fermi energy in the thermodynamic 
limitQQ. Though this minor difference in definition does not affect the results and the physics, Tk defined here offers 
an intuitive understanding of the crossover at finite sizes corresponding to case (I) and (II) mentioned above-for 
K < Kc, Tk is expected to increase with increasing size L for case (I) as the system recovers the Kondo resonance in 
Loo limit, and for K > Kc Tk is expected to decrease to as the system approaches to the local spin-singlet ground 
state in the thermodynamic limit (see details below). Nevertheless, we have checked that with the definition in Ref. 
and our results for Tk at finite sizes indeed reproduce those in Ref.Q and @. We have also checked that our 
results on Tk for K = are consistent with the finite-size behaviors in Egg and E^^ via perturbation theory in Ref. 

a. 

In the Kondo phase, Tk reduces to the Kondo temperature for the single quantum dot embedded in a mesoscopic 
ring {K — 0). When L oo,Tk in the Kondo phase approaches to T^. In the crossover between Kondo and the local 
singlet phases, Tj. decreases with increasing K , and finally — > when the system is in the local singlet ground state 
where the Kondo effect is completely suppressed by the antiferromagnetic spin-spin couping. We analysize below in 
details the crossover between the Kondo and the local spin-singlet ground states from the behaviors of the Kondo 
temperature. 



1. Varies K with fixed L 
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FIG. 3: Tk as functions of K. (a)Case I: L = 4n (b)Case 11: L = 4n -I- 2. Other parameters: = -0.8t, Ir = II = 0.4t, = 0. 
Here we set t = 1 as the unit. 

We first monitor how the Kondo resonance is destroyed by varying the antiferromagnetic spin-spin coupling strength 
K &t a. fixed finite size L. Fig.Q shows Tk as a function of antiferromagnetic coupling strength if at a fixed size 
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L = 300 (or ~ 0.3). In both case (I) and (II), Tk vanishes as iiT ^ Kc2, indicating the suppression of Kondo 

resonance by the antiferromagnetic spin-spin interaction. However, there are minor differences between these two 
cases in how fast Tf. vanishes as K close to Kc2 ■ In case (II) , T^. remains a constant over a wider range of i^T < Kd 
compared to that in case (I) before it decays. This suggests that the Kondo effect at a finite size seems more robust 
in case (II) (L = 4?! + 2) than in case (I) [L — An, An + l,4n + 3) so that the crossover region in case (II) is much 
narrower than in that in case (I). This is consistent with our numerics as Tk for i = 4n + 2 is found to have the 
largest value among L (mod 4) for a given size. 



2. Varies L with fixed K 



Next, we present the results on the finite-size dependence of the Kondo resonance at a fixed antiferromagnetic 
spin-spin coupling strength K. In analogous to the Numerical Renormalization Group (NRG) method, the ground 
state is computed and monitored as we decrease the energy scale 1/L (or equivalently increase the system size L) 
until we reach the thermodynamic limit L — > oo (or effective zero temperature). 

First, we describe the qualitative behaviors of Tk at finite sizes in the two cases as mentioned above. For case (I) 
(represented by i = 4n, see Fig|4Ka)) Tk increases with increasing L as the crossover region is for K < Kc where the 
Kondo resonance is recovered at large system size L ^ ^j?; while for case (II) {L = An + 2, see FigHJb)) since the 
crossover occurs for K > Kc, Tk decreases with increasing L and it vanishes in the thermodynamic limit. It is clear 
from the crossover behavior in Fig. |4]that the finite-size effect appears for Ll^^ where the Kondo screening cloud has 
not yet fully developed for K < Kc and has not yet completely destroyed for K > Kc', while this effect diminishes as 
the system approaches to the thermodynamic limit or L Note that in case (I) with large K > Kc and case (II) 

with small K < Kc, the ground state remains at local singlet and Kondo state, respectively; therefore, no crossover 
behaviors are found. It is also worthwhile noting that Tk for L = 4n -I- 2 (case (II)) ioi K = approaches to T° 
from above as L approaches to the thermodynamic limit. This suggests that the Kondo effect in this case seems more 
robust at finite size than in the thermodynamic limit. Therefore, at any finite size L, it is necessary to apply a larger 
antiferromagnetic coupling K compared to Kc which is required in L — > cxd limit to suppress the Kondo effect. This 
provides an explanation why we always find the crossover behavior for i = 4n -t- 2 for K > Kc- On the other hand in 
case (I) {L — An, An -I- 1, 4n + 3) for K = 0, Tk approaches to T^ from below as L increases, which explains why the 
crossover occurs for K < Kc in this case. 
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FIG. 4: Tk as functions of 1/L: (a) L = An, K varies from 0.0025f to 0.015f, distance: 0.0025t. (b) L ^ An + 2, K varies from 
0.036t to O.OAt, distance: O.OOlf. Other parameters: ed — -0.8t, tn — II = 0.4t, = 0. Here, we set t = 1 as the unit. 



To investigate the nature of the quantum phase transition in the thermodynamic limit more closely, we then perform 
the finite-size scaling for Tk in the crossover region. We find Tk in case (I) and (II) follows its own unique universal 



8 



scaling function of 1/{T*L) as K ~> Kc (see Fig. [5]), where T* is the crossover energy scale defined as: 

T* ^cfkei,p{~7rfk/ \K -K,\) (13) 

with Tk = c'Tk has the same form in both case (I) and (II). Here, c, c', are non-universal fitting prefactors 
depending on the initial parameters of the Hamiltonian. For ed = -0.8t, tji — — 0.4t, </) = 0, wc find c « 5.5, 
Kc « 0.0271 (in unit of t) in both case (I) and (II); c' w 0.65 in case (I), and c' « 0.5 for case (II). Note that the 
expression for T* in Eq. [13] is quite general as K can be either smaller (case (I)) or larger (case (II)) than K^- The 
universal scaling at finite sizes and the same exponential form for the crossover energy scale T* valid for both cases 
strongly indicate that in the thermodynamic limit the system exhibits the Kosterlitz-Thouless transition at a finite 
critical antiferromagnetic spin-spin coupling strength > 0. We have checked the consistency of our result from 
our finite-size scaling that Kc indeed reaches to the same value in the thermodynamic limit for both case (I) and (II) 
even though the corresponding crossover regions are on the opposite side of the transition {K < Kc for case (I) and 
K > Kc for case (II)). 

Note that unlike the similar side-coupled double-dot system studied previously in Ref . [3] where the KT transition 
between the Kondo and spin-singlet phase occurs at Kc = 0, we find in our setup a finite Kc > for the same KT 
transition. We think this difference might be related to the more singular DOS of the conduction electron bath in the 
current setup of a tight-binding ring compared to that in a continuous Fermi sea with a constant DOS, making the 
Kondo resonance more robust against the antiferromagnetic spin-spin coupling in the current set up and therefore 
leads to a finite Kc instead of Kc — 0. However, it is also possible that the finite Kc we find here is the artifact of 
the large-N slave-boson mean-field theory. Further study is necessary to clarify this issue. 

Finally, as a consistency check, as the system approaches to the thermodynamic limit, L —> oo, Tk ^ for L = An 
(case (Vj, K < Kc); while ^ for L = 4n -f 2 (case (II), K > Kc), as expected (see Fig. [5|). 
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FIG. 5: Tk as functions of l/(r*L): (a) L ^ 4n, K varies from 0.0025* to 0.015t, distance: 0.0025*. (b) L = 4n + 2, K varies 
from 0.036t to 0.04t, distance: O.OOlt. Tk follows universal scaling functions of 1/{T*L). Other parameters: ed = -0.8t, ta = 
tL = 0.4t, <j!> = 0. Here, we set t — 1. 



C. Local Density of State 

We now turn our attention to the local density of states (LDOS) on the dot 1, given by 

Pqd{co) ^--ImG^d^dA^) (14) 

TT 

where ^ (t) =< di{t)d\{0) > is the retarded Green's function of dot 1 which directly couples to the ring. Near the 
Fermi surface, pqd{lj) determines the transport properties of the system, and it can be obtained from the Green's 
function Gdidi (f-^) via equation of motion approach. First, the mean-field Hamiltonian in momentum space is given 



9 



by: 

Hmf = ^emcLc,„ + i/.C. + e~d/t^/i,-^(i„4/i+i7.C.) (15) 

m,a m 

- x/l/2 - H.C. + ^ + Xibl - 1) 



Where 



/2 m ~ ■ 
-sin(-^)[tL+i^exp(#)(-l)"+i] (16) 

with m = 1,2, • • - L - 1^. 

The equation for a general retarded Green's function G^{oj) is then given by: 

{lo + ir^ - HMF)Gf jiio) ^ I (17) 

where / is the identity matrix, i,j can be /i, /2, or m. We therefore get the following three equations for Gf^f^, Gf^m 
and G/,/^: 

(c + zr, - e~rf)G«/,(c.) - xGf,y, -Y^t^G^^^ioj) = 1 (18) 



-xGf,,, + (^ + z^)Gf,,,=0 (19) 

-t„^Glf^+{u + ii^~e^)Glf^^Q (20) 
The above equations are easily solved and we get: 

G^ , iuo) =hlG^ f = ° , , „ (21) 



ed + - E 



where r/ 0. We then obtain pqd{llj) by Eq. [TH In the following we analysize the crossover behaviors in LDOS 

1. K = Q 

In the absence of the antiferromagnetic spin-spin coupling {K = 0) the LDOS has been extensively studied where 
Pqd{oj) depends sensitively on L (mod4) at finite sizes @. As shown in Fig.®, our results on LDOS in this case 
at finite sizes are qualitatively in good agreement with that in Ref.f^. For the convenience of later discussions we 
summarize below the behaviors of LDOS in the four different cases of L. The key features are- (i). There exsits a 
main Kondo resonance peak (located either symmetrically or asymmetrically with respect to cj = 0) followed by pairs 
of side peaks, (ii). As shown in Fig. [SI as L increases the main Kondo peak gets more pronounced and closer to 
w = 0; while the side peaks are gradually merged into the broadened main peak. In particular, the LDOS for L = An 
is very symmetric pqd{(^) — Pqb{~^\ suggesting the symmetry between particle and hole excitation energy in the 
finite size spectrum. The asymmetric Kondo peaks for L = An + I and i = 4n + 3 are located on the opposite sides 
(left for L = 4n + 1 and right for L = 4n + 3) of uj. For L = An + 2, however, the LDOS shows asymmetric double 
Kondo peaks with comparable sizes and a dip at = 0. These differences in DOS among the four sizes of L (mod 4) 
can be explained in terms of the energy levels corresponding to the highest occupied (HO) and the lowest unoccupied 
(LU) states relative to the Fermi level[5|]: For L = An, both HO and LU levels are around the Fermi energy ep = 0, 
leading to a symmetric single peak in LDOS at u; — 0. For L = 4n ± 1, HO and LU levels are both below and above 
ep, respectively, giving rise to an asymmetric single peak in LDOS below and above uj = separately. However, for 
L = 4n + 2, LO and HU levels are on the opposite side of the Fermi level, resulting in splitted double Kondo peaks 
below and above ui = 0. It should be noted that despite the differences at finite sizes {L ~ ^°), the LDOS in all four 
cases indeed merges to a single Kondo peak as the system reaches the thermodynamic limit L ^ oo (or L ^ ^°). (iii). 
The LDOS on dot 1 obeys the following symmetries [5| : Pqjj{lo, 4>) = p]^{uj, (/) + 7r) and pq]j{~uj, 0) = pqd{^-, 4> + it). 
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FIG. 6: Local density of state for K = system (in arbitrary unit) (a) L — 4n (b) L — 4n + 1 (c) L — An + 2 (d) L — An + 3. 
Here, r — £,^./L, the ratio between and L. Other parameters: = -0.8t, tu — II = 0.4t, (j) — Q. Here, we set t = 1. 

We would like to point out here that from the width D of the central Kondo peak(s), which can be approximately 
regarded as a quantity proportional to the effective Kondo temperature T^, in LDOS ior K = Q case one can 
qualitatively understand the opposite trends in Tk at finite sizes in case (I) and (II) mentioned above. For L = 
An,4,n + l,4n + 3 (case (I)), the width D becomes larger as L increases, indicating an increase in as the system 
approaches to the thermodynamic limit; while as ior L — 4n + 2 (case (II)), D gets smaller as L increases, suggesting 
a decrease in Tk as L ^ oo. We analysize in details below the behaviors of LDOS for K > 0. 

2. K >0 

At a finite K > 0, the LDOS on the dot shows a crossover between the Kondo phase and the spin-singlet phase. For 
K < Kci at a fixed size L, the LDOS remains the same as that for K = Q. For L = 4n, 4n + 2 we find a continuous 
evolution in LDOS from Kondo to the crossover region near K = Kd] while for L = 4n + 1 , 4n + 3 the LDOS exhibits 
a first-order jump at Kd- For K > Kc2, we find pqd{ljj) = as the indicator of the local spin-singlet phase since 
&o = 0. In the crossover region Kd < K < Kc2, the Kondo peak in LDOS splits into two with respective to w = 0. 
The splitting gets wider as K increases further. Details are shown below. 
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(i) . Varies K with fixed L: 
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FIG. 7: Local density of state (in arbitrary unit) with various Kiorr^ ^l/L ~ 0.53. (a) L = 4n (b) L = 4n + 1 (c) L : 
(d) L = 4n + 3. Other parameters: td = -0.8t, tu = II ~ 0.4t, = 0. Here, K is in unit of t, and we set t = 1. 



4n + 2 



In FiglJl we show LDOS for various antiferromagnetic spin-spin coupling strength if at a fixed size L w 200 
(5°/L ^ 0.5). The behaviors of LDOS are classified by L(niod 4). The common features in Fig. [7] are as follows. 
In each of the four LDOS plots in Fig. [71 the Kondo peak in pq]j{llj) splits into two at a small value oi K > K^i, 
indicating the crossover between Kondo and local spin singlet phases. The two peaks separate further apart and 
become more symmetric as K increases. At the end, pqd{^) vanishes when K > Kc2 in the local spin singlet phase. 
Note that the values for Kc2 depend sensitively on L (mod 4). This can be understood as when L « 200 the effective 
Kondo temperature Tk in case (II) is much larger than that in case (I) (see FigHl Fig. [5|) , which explains why we 
need a larger value of Kc2 in case (II) to suppress the Kondo effect than that in case (I). We present details below. 



(ii)- Varies L with fixed K: 



Fig. ([8]) shows how pq£i{lo) changes with the system size L. For < K < Kc, starting from a small size L < Ld, 
the LDOS changes from the behavior of a local spin-singlet state with pqd{uj) = to that in the crossover region 
with splitted peaks and finite LDOS at w = 0, and finally it recovers the Kondo resonance at larger size L > Lc2- 
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FIG. 8: Local density of state (in arbitrary unit) with various L and fixed K (in unit of t). Here, r = ^°/L. (a) L = 4n (b) 
L = 4n + 1 (c) L = 4n + 2 (d) L = 4n + 3. Other parameters: = -0.8t, tn = II ^ 0.4t, (f) = 0. Here, we set t ^ 1. 



We separate the discussion here into two different cases: K < Kc (case (I)) and K > Kc (case (II)). 
Fig.(I5I^a),(b),(d))show finite size dependence of pqd{^^) &t K < Kc (case (I)) for L — 4n, An + 1, and An + 3, 
respectively. For small system sizes, the system is in the local spin-singlet state with vanishing LDOS. At intermedi- 
ate sizes in the crossover region, however, LDOS starts to develop peaks with a dip at w = 0. As size further increases, 
these splitted peaks either gradually {L = An) or suddenly (L = 4n + 1 and An + 3) merge into a single Kondo peak 
located either symmetrically (L = An) or asymmetrically {L — An + 1 and An + 3) with respect to w = 0. These 
Kondo peaks then follow the evolution at finite sizes for K — and finally recover the single symmetric Kondo peak 
in the thermodynamic limit. 

For K > Kc (case (II)), the size dependence of pqaii^) is shown in Fig.ilHI^c)) with K — 0.03i > Kc- The LDOS 
exhibits a crossover from the Kondo to local spin-singlet state with increasing size L. For K ~ Kc the LDOS changes 
from singlet state at small sizes to the crossover regime within the range 100 < L < 800 that we investigate. To study 
the the ultimate fate of the system one needs to go much larger system size than L ^ 800, which is beyond the scope 
of our computational limit. Note that our analysis on the LDOS by changing K with fixed L (changing L with fixed 
K) is consistent with the green line (red line) of the schematic phase diagram shown in Fig. ([2]). Finally, the first 
order jumps seen in LDOS a,t K = Kd for L = An + 1 and L ^ An + 3 may be due to the artifact of the mean-field 
theory. Further investigation is needed to clarify this issue. 



D. Persistent Current 



We now analysize the crossover from the behaviors in persistent current. The Kondo screening cloud in our closed 
setup is restricted itself in the ring. To get an experimental access, one possible way is to measure persistent current 
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(PC) induced by changing the magnetic flux threading the ring without attaching leads to it. PC is defined as: 



h del) 



(22) 



PC can be served as a detector for the Kondo screening cloud since as a function of (j) PC behaves very differently 
when ^fc ^ L than when ^fc <g; L. Such persistent current experiments have been reported recently on micron sized 
ring without containing the dot[l3|- In our setup, PC can be used to measure the strength of the Kondo correlation 
as the antiferromagnetic coupling K is expected to suppress the Kondo effect, leading to a weaker PC. 



K=0 




FIG. 9: Persistent current (in unit of Jo) versus magnetic fiux (in unit of 27r/$o) with K = Q. Red dashed line is for 
r = ^° ~ 0.53. Black solid line is for ^l/L ~ 1.06. (a) L = 4n (b) L = 4n + 1 (c) L = 4n + 2 (d) L = 4n + 3. Other parameters: 
Ed = -0.8t, tn — tL = 0.4t. Here, we set t = 1. 

Before we present our new results on the PC, it proves useful to summarize the previous results[3]@| for K = 0. 
As shown in Fig.®, since the system is in the Kondo regime, PC increases with increasing system size L. PC as a 
function of the magnetic flux is a weak sinusoidal wave when L <^ while it behaves like a saw-tooth when L ^ 
Meanwhile, there exists a relation between PC for L and L + 2: In{4>) ~ -^Af+2(0 + 7r)[3] |1| due to the tt shift in LDOS 
of the quantum dot from the system with N sites to iV + 2 sites. Furthermore, the magnitude of PC for L = odd is 
much smaller than that for L = even. This can be understood as for L = odd the electrons are fully occupied at the 
Fermi level, which suppresses the Kondo resonance; while as for L — even there is an unpaired electron on the Fermi 
level, leading to the Kondo resonance. The results ioi K > are shown below. 
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1. Varies K with fixed L: 




FIG. 10: PC versus (j) with various for L ~ 240 {^l/L ~ 0.44). K = 0,0.015, 0.03 (in unit of t) (a) L 
(c) L = 4n + 2 (d) L = 4n + 3. Other parameters of data: td = -0.8t, Ir = II = 0.4t. Here, we set t = 1. 



4n (b) L = 4n + 1 



PC (in unit of Jq where /q = evp/L is the persistent current of an ideal metalhc ring with vp being Fermi velocity) 
versus magnetic flux with different antiferromagnetic coupling strength i^T at a finite size L ~ 240 {^^/L ~ 0.44) is 



shown in Fig. (|10p . The general trend in all four cases is that the amplitude of PC gets smaller as K is increased 
and finally vanishes for K > Kc2, in agreement with the expectation on the suppression of the Kondo effect by the 
antiferromagnetic spin-spin coupling. Also, as shown in Fig. llOi the shape of the PC as a function of 4> changes from 
saw-tooth shape at smaller K values to weak sinusoidal waves at larger K values, similar to that for K = 0. There are 
detail differences among the four cases. In Fig. (|10r a) and (c)), with increasing K the persistent current /l=240 near 
O.Stt < (p < I.Stt decreases more slowly than that for On < cj) < O.Stt and I.Stt < (p < 27r; similarily for /l=242 except 
for the range of cj) is being interchanged. This suggests that the Kondo effect is more robust in O.Stt < <j) < I.Stt for 
L — 4n and in On < cj) < 0.5n and I.Stt < (j> < 2n for L = 4n + 2. Note that we find the similar jumps at K = Kci 
mentioned previously in LDOS to also appear in PC for L = An + 1 and L — An + 3; while PC is continuous at at 
K = Kci for the other two cases. Despite the above detail differences, we find two common features in Fig. (jTU]) that 
remain the same as in the case oi K — 0: (i). In{4') — In+2{4' + "') a-nd (ii). PC for L = even is larger than that for 
L = odd. 
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2. Varies L with fixed K: 



Fig. fTTj) and Fig. (fT2|) show PC versus magnetic flux at different system sizes. In the foUowing the discussion is 
separated into two cases: (i). K < (case (I)) where PC increases with increasing size L, and (ii). K > (case 
(II)) where PC decreases to as L reaches the thermodynamic hmit. 

(i). K<K, 
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FIG. 11: Persistent current with various L and fixed K (in unit of t). {a.) L = An + 1, K = 0.008 (h) L = An + 3, K = 0.008 
(c) L = 4n + l, K ^ 0.015 (d) L ^ 4n + 3, K ^ 0.015. Here, r = ^^/L. The parameters are: e<j = -0.8t, Ir ^ II ^ 0.4t, = 0. 
Here, we set t = 1. 

We first look at cases with K < Kc- Fig. pTj) shows PC versus magnetic flux with different size L for K = 0.008 
and K = 0.015 (in unit of t). For K = 0.008 (see Fig[TTJa)(b)), as L is increased we find the PC exhibits two different 
behaviors: near (f> ~ n for L ~ An + 1 and = 0, 27r for L = 4n + 3 PC changes from the behavior in crossover region 
to the Kondo state at if = (see Fig. (|Hb)(d)); while out of these ranges of (j) its behavior remains the same as in 
the crossover region but with an increasing amplitude. A first-order jump is seen to separate these two regions. We 
expect at much larger system size the Kondo phase will be restored eventually over the entire range of cf). We have 
also investigated PC for L = An, An + 2 (not shown here) and find the same qualitative behaviors except that instead 
of jumps we find a continuous change in PC. Also, the behavior in I^n+i (or equivalently if^+s) is the same as that 
for if = for all the sizes we investigate, suggesting the Kondo phase is very easily restored in these ranges of (f>. 

On the other hand, for a larger value of if = 0.015 (see FigfTTTcHd)). as L increases from L ~ 60 to L 420 
(or r = £,1/L changes from 1.77 to 0.25) PC still stays in the crossover region with an increasing amplitude. This is 
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expected as for larger value of K one must go to much larger system size to observe the restoring of the Kondo effect 
in PC. 

(ii). K>K, 
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FIG. 12: Persistent current with various L and fixed K. 
data: = -0.8t, Ir = II = 0.4t, K = 0.036 (in unit of t) 



(a) L = 4n (b) L = 
Here, we set t = 1. 



An + 2. Here, r = Other parameters of 



We now investigate the case where K > K^- As shown in Fig. p^ . since K > in this case, PC decreases in 
amplitude with increasing system size L and finally the system reaches the local spin-singlet state with vanishing PC. 
Note that if!^^ and correspondingly lf^^2 ^^'^ always vanishingly small for all the sizes we investigate, suggesting the 
systems are already in the local spin-singlet states from the start for these ranges of (j>, which is consistent with our 
mean-field phase diagram. However, for the remaining ranges of (j) the systems start from the crossover region for 
smaller sizes and reach finally to the local spin-singlet state at large size. 

From above results for PC, it implies that by applying magnetic flux, we may change the ground state of our system 
at finite sizes either from the local spin-singlet state to the crossover region or from the Kondo phase to the crossover 
region. 



IV. Conclusions 



To summarize, we have studied via large-N slave-boson mean-field approach the Kondo effect in a side-coupled 
double-quantum-dot system where one dot is embedded in a mesoscopic ring. The competition between the Kondo 
effect and the antiferromagnetic spin-spin interaction in this geometry gives rise to the Kosterlitz-Thouless quantum 
phase transition with a finite critical value Kc in the thermodynamic limit. The mean-field phase diagrams of the 
model depends on the finite size L (mod 4): for L = 4n, 4n-|- 1, 4n-|-3 (case (I)) the crossover occurs for K < Kc, while 
for L — An + 2 (case (II)) the crossover exists for K > Kc- To further study how the Kondo screening is suppressed 
by the RKKY, we have performed a systematic finite-size analysis on the Kondo temperature Tk, the local density 
of states on the dot pqd{uj) which connects to the ring , and the persistent current PC induced by the magnetic flux 
penetrating through the ring. For a fixed K < Kc, we find all the above quantities flow to the Kondo phase (same 
phase as in JsT = 0) with increasing system size: Tfc and PC increase in magnitude and pqd{u) develops a pronounced 
single Kondo peak centered at = 0; while for K > Kc, the flow at finite sizes is towards to the local spin-singlet 
phase where all of the three observables decrease in magnitude to with increasing size. From the finite-size scaling 
of Tk, we have shown that Tk is an universal function of 1/{LT*) for 1/L < T* where T* has been unambiguously 
identified the characteristic Kosterlitz-Thouless crossover energy scale: T* = cTkexp[—TTTk/{K — Kc)]- For a fixed 
size L, we find the ground state remains in the Kondo phase for K < Kd and it crossovers to the local spin-singlet 
phase for Kd < K < Kc2 and finally reaches the spin-singlet phase for K > Kc2- For L — An + l,An + 3, we find 
the first order jumps in all of the three observables at the phase boundary between Kondo and the crossover region. 
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Whether these first order jumps are the artifacts of the large-N mean-field theory is yet to be clarified. Note that 
unlike the similar system studied previously on the side-coupled double-dot system embedded in conduction electron 
Fermi sea with constant density of states, the key finding in this paper is the KT transition with a finite critical 
value Kc > in a side-coupled double-dot system embedded in a mesoscopic ring (in stead of K^, = for the system 
with Fermi sea of continuous spectrum). Whether this finite Kc is due to the artifact of the mean- field theory or 
due to the more singular density of states of the ID tight-binding ring needs further investigations. Our results on 
the transport properties of the system are relevant for future experiments on the side-coupled double-quantum-dot 
system embedded in a mesoscopic ring. 
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